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We study the time evolution of correlation functions, spin current, and local magnetization in an 
isolated spin- 1/2 chain initially prepared in a sharp domain wall state. The results are compared 
with the level of spatial derealization of the eigenstates of the system which is measured using 
the inverse participation ratio. Both integrable and non-integrable regimes are considered. Non- 
integrability is introduced to the integrable Hamiltonian with nearest neighbor couplings by adding a 
single site impurity field or by adding next-nearest-neighbor couplings. A monotonic correspondence 
between the enhancement of the level of derealization, spin current and magnetization dynamics 
occurs in the integrable domain. This correspondence is however lost for chaotic models with weak 
Ising interactions. 
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I. INTRODUCTION 

The dynamical behavior of quantum many-body sys- 
tems out of equilibrium offers several open questions. In 
terms of transport properties, for example, the conditions 
for ballistic or diffusive transport have received a good 
deal of attention throughout the years [TI-fTll] . The inter- 
est in the subject has been further reinforced by exper- 
iments with magnetic compounds, where ballistic trans- 
port of heat has been reported [73,1171. Another stimulat- 
ing scenario appears in the context of experiments with 
optical lattices 0, [l5| , where the dynamics of strongly 
correlated quantum systems can be investigated for long 
times. Here, the possibility of studying the pro blem of 
thermalization in isolated quantum systems [16j has in 
part motivated the renewed interest in this topic of re- 
search mm. 

In the present work we are interested in the quench dy- 
namics of one-dimensional systems of interacting spins- 
1/2. We focus on a particular initial state corresponding 
to a sharp domain wall where the spins point up in the 
first half of the chain and down in the other half. The 
experimental realization of such inhomogeneous states re- 
quires a magnetic field gradient, which has been realized 
in an optical lattice loaded with atoms in two hyperfine 
states [19| ■ The domain wall state was considered before 
in numerical studies of magnetization transport in anti- 
ferromagnetic spin- 1/2 XXZ chains J2(j, where using the 
adaptive time-dependent density-matrix renormalization 
group method, a sharp suppression of transport was ver- 
ified for large interactions. It also appeared in studies of 
the dynamics of a spinon in an external magnetic field [2l[ 
and was used in numerical comparisons of transport be- 
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havior in integrable and nonintegrable systems [22M24| . 
In the context of thermalization, the dynamics of a sharp 
domain wall in the gapped phase of the XXZ chain was 
studied using Algebraic Bethe Ansatz in [25]. 

The evolution of a domain wall state was investigated 
analytically in [26| and [27j for the XX model and numer- 
ically employing the truncated Wiener approximation for 
the quantum sine-Gordon model [281 ] . The domain wall 
in the XX model and the Luttinger liquid was found 
to spread out ballistically into a current carrying steady 
state [I?], . In a spin-chain, current flow implies a ro- 
tation in space of the transverse component of the spins. 
As a result, the transverse spin correlation function in 
these simple models were found to reach a nonequilib- 
rium steady state which showed spatial oscillations at 
a wavelength inversely related to the magnitude of the 
current [27|. As back-scattering interactions were added, 
the current was found to decrease, and consequently, the 
wavelength of the oscillations in the transverse spin cor- 
relation function was found to increase (28|. These re- 
sults are further explored here taking into account both 
integrable and nonintegrable systems. The approach is 
to study the time-evolution numerically employing EX- 

POKIT min. 

EXPOKIT is a software package based on Krylov sub- 
space projection methods. Instead of diagonalizing the 
complete system Hamiltonian H, we use this package to 
directly compute the action of the matrix exponential 
e -tHt on a vec t or f interest. We consider spin-chains 
with short range couplings and 20 sites leading to sparse 
matrix Hamiltonians of dimension 184 756. 

Based on transverse spin correlation function, spin 
current, and local magnetization, we perform a numeri- 
cal comparison between integrable spin- 1/2 chains with 
nearest-neighbor couplings only and their nonintegrable 
counterparts. The latter are realized by adding a single- 
site impurity field or by adding next-nearest-neighbor 
couplings. The wavelength of the oscillations of the trans- 
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verse spin correlation function increases with the Ising 
interaction for the three models. More unexpected are 
the results for current and magnetization for the chaotic 
models. For the integrable model, the behavior is pre- 
dictable. As the Ising interaction increases, the eigen- 
states of the system become spatially more localized, 
which reflects in the motion of the excitations: the spin 
current decreases and the decay of the absolute local 
magnetization slows down. For large interactions, the 
same tendency is verified for the chaotic models, whereas 
surprises appear as the strength A of the Ising interac- 
tion changes from zero to small values. For the model 
with impurity, the addition of a small interaction leads 
to the onset of chaos, which is followed by a significant 
increase of the level of spatial derealization of the eigen- 
states. Contrary to that, the magnitude of the current 
and the decay rate of the magnetization decrease. In 
the case of next-nearest-neighbor couplings, the system 
is chaotic even in the absence of Ising interaction (pro- 
vided that the ratio between next-nearest- and nearest- 
neighbor couplings is sufficiently large). Adding A con- 
tracts the cigenstatcs, whereas, quite countcrintuitively, 
the spin current increases and the decay of the absolute 
local magnetization is enhanced. 

For the model with impurity, the opposite behavior 
between derealization and dynamical spreading of exci- 
tations may be an indication of the transition from bal- 
listic to diffusive transport. We do not have a complete 
explanation for the behavior of the chaotic system with 
next-nearest-neighbor couplings, but hope this work will 
motivate further investigation. 

The paper is organized as follows. In Sec. [Til the mod- 
els and the physical quantities studied are defined and 
a description of the numerical method is provided. In 
Sec. IIII1 the results for the correlation function, magne- 
tization and current are analyzed, and in Sec. IIV1 these 
results are compared with the level of spatial derealiza- 
tion of the eigenstates of the system which is measured 
in terms of the inverse participation ratio. Concluding 
remarks are presented in Sec. [V] 



II. MODELS AND NUMERICAL METHOD 

We study a one-dimensional spin- 1/2 system with open 
boundary conditions described by the Hamiltonian 



H — H z + i?NN + O^NNN, (1) 

where 



H z - e L/2+lSl/ 2+ i, (2) 
L-l 

#NN = ^ J + ^n^n+l) + ^n^n+l] > (3) 

n=l 
L-2 

#NNN = 2_j J [( S n S n+2 + S n S n+2) + AS£/S£ +2 ] ■ ( 4 ) 
n=l 

Above, h is set equal to 1, L is the number of sites, and 
g%,v,z = a^' v,z /2 are the spin operators at site n, u^' y ' z 
being the Pauli matrices. 

(i) The system is clean when all Zecman splittings, 
determined by a static magnetic field in the z direction, 
are equal and it has on-site disorder when the spins have 
different energy splittings. We will consider the case of a 
single-site impurity field where only the Zeeman splitting 
e L/2+i a t site L/2 + 1 is different from zero. 

(ii) The XXZ Hamiltonian i?NN describes the nearest- 
neighbor (NN) exchange; S*S* +1 is the Ising interaction 
and S%S% +1 + S^S^ +1 is the flip-flop term which tends to 
orient the spins in the XY plane. The coupling strength 
J and the anisotropy A are assumed positive, thus favor- 
ing antiferromagnetic order. By varying A in H = -Haw, 
the system moves from the gapless XY phase (A < 1) to 
the gapped Ising phase (A > 1), which, at T = 0, cor- 
responds to the metal- insulator Mott transition [3l| . In 
the particular case of A — 0, we deal with the XX model; 
whereas A = 1 corresponds to the isotropic Heisenberg 
Hamiltonian (32[. 

(iii) The parameter a refers to the ratio between the 
next-nearest-neighbor (NNN) exchange, as determined 
by -Hnnn i and the NN couplings. The inclusion of NNN 
antiferromagnetic exchange frustrates the chain, since 
NN exchange favors ferromagnetic alignment between the 
second neighbors. 

Depending on the parameters of Hamiltonian (JT|), dif- 
ferent symmetries are found. The total spin in the z 
direction, S z = X)n=i ^ s conserved for all parame- 
ters. When cl/2+i = 0, other symmetries include [33| : 
invariance under lattice reflection, which leads to parity 
conservation; invariance under a 7r-rotation around the 
x-axis in the case of S z = 0; and conservation of total 
spin, S 2 = (ELi Sn) 2 when A = 1. 

A. Integrable and chaotic regimes 

Adjusting the parameters of the Hamiltonian, the sys- 
tem may undergo a crossover to the chaotic regime. We 
consider here the following scenarios. 

• Integrable model with cl/2+i = 0.0, a = 0.0; H = 
-ff/vjv- A clean one-dimensional spin-1/2 model with 
NN couplings only is integrable. When A = 0, it 
can be mapped onto a system of noninteracting spinlcss 
fermions [3J]. When A ^ 0, its ground state properties 
can be studied using Bethe Ansatz (3f|. 



3 



• Chaotic model with £l/2+i 7^ 0, a = 0.0; H = H z + 
Hjfjf. The addition of on-site disorder (3(| to the XXZ 
model, even if it corresponds to a single impurity in the 
middle of the chain [37|, leads to the onset of quantum 
chaos provided £l/2+ i ^ J (if the impurity becomes too 
large it breaks the chain in two pieces) and A < 1, but 
not zero. The presence of interaction is essential; it is 
the interplay between interaction and disorder that leads 
to chaos. Here, we will consider an impurity of strength 

(-L/2+1 = 0-5 J. 

• Chaotic model with £l/2+i = 0.0, a 7^ 0; H = 
Hnn+oiHnnn. The inclusion of NNN couplings may also 
drive the system to the chaotic domain for sufficiently 
large values of a f38M40| . We will consider a = 0.4, 
which guarantees chaoticity when A < 1 . It is interesting 
that the system is chaotic even when A = 0, which we 
verified by computing the level spacing distribution [4l| . 
When a 7^ and A = 0, the Jordan- Wigner transforma- 
tion does not map the spin- 1/2 chain onto a system of 
nonintcracting spinless fermions, but instead correlated 
hopping terms appear, which is the cause for the onset 
of chaos. The level spacing distribution is the most com- 
mon way to verify whether a system becomes chaotic. It 
requires the complete diagonalization of the Hamiltonian 
and, very importantly, the separation of the eigenvalues 
according to symmetry sectors. If eigenvalues from dif- 
ferent subspaces are mixed, we may obtain a Poisson dis- 
tribution (indicating integrability) even if the system is 
chaotic. 



B. Quantities of interest 

In Sec. IIII1 we will study the dynamical properties of 
three quantities of interest, which are described below. 

(i) The transverse spin correlation function corre- 
sponds to 

C xx {n,n + m) = {S x n S x n+m ). (5) 

(ii) The local spin current, I s ,m is associated with 
the conservation of total spin in the z direction, S z = 
Sn=i Sn an d obeys the continuity equation, 

—^ + div(I s , n )=0. 
at 

In the bulk [Hj], 

-i[H, S z n ] = div{I s , n ) = (/,,„ - I s ,n-i), 
but in the extremes, since the chain is open, we have 
-i[H,S*] = div(I s j) = I sJ 

and 

-i[H,Sl] = div(I a , L ) = -I a , L - l . 



From the equations above, we find that the local spin 
current for a system with NN couplings only, with or 
without impurity, is given by 

where 1 < n < L — 1. 

In the case of the model with NN and NNN couplings, 
the local spin current becomes 

for 2 < n < L — 2, and at the borders, 

= J{SIS V 2 - SfSZ) + a,J(S^Sl - SfSI), 

Equations. (J5J) and (J5]) imply that for the XXZ chain 
a nonzero current gives rise to a static spin configuration 
where the spins rotate in the XY plane from one site to 
the next at a rate which increases with the magnitude 
of the current. However, as is apparent in Eq. ([7]), this 
simple relation between the NN spin correlations and the 
current is lost for NNN couplings. 

(iii) The local magnetization of each site is, 

M„ = <*(«|*(*)>, (8) 

while the local magnetization of the first half of the chain 
is given by, 

L/2 

M {1>L/2) = (*(t)\j2s z n m))- (9) 

n=l 

The results for the dynamical behavior of the above 
quantities will then be compared with the level of der- 
ealization of the eigenstates in Sec. IIVI 

C. Numerical method 

We use EXPOKIT [H, HI to study numerically the 
time evolution of the system. This is a software pack- 
age that can deal with large sparse matrices, as the ones 
we have in this work. It avoids the complete diagonal- 
ization of the Hamiltonian and instead provides routines 
to compute the matrix exponential e~ lHt . The algorithm 
computes the matrix exponential times a vector 1^). The 
main idea is to approximate e - ^'!^) by an clement of 
the Krylov subspacc of dimension k, which is small com- 
pared to the dimension of the total Hamiltonian. The 
Krylov technique is based on the Arnoldi iteration when 
dealing with a general matrix and on the Lanczos it- 
eration if the matrix is symmetric or Hermitian. The 
product e~ lHt \^>) is approximately evaluated within the 
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low-dimensional Krylov subspace. The problem is then 
reduced to the calculation of an exponential of a dense 
k x k matrix, which is done by using the rational Pade 
approximations. 

Here, we consider k = 30 and time steps of 0.01/ J. The 
software proved to be very fast and accurate. For a chain 
with L = 16, for example, we computed the deviation 
\M exact — M expo kit \ between the result for site magneti- 
zation obtained from exact diagonalization (M exact ), and 
that obtained with EXPOKIT (M expo kit), for all sites of 
the chain. This was done for all 3001 instants of time be- 
tween and 30/ J. For a system with 8 up-spins, A = 0.5, 
and for both NN and NNN couplings (a = 0.4), the de- 
viation was always smaller than 10~ 8 . Using complete 
exact diagonalization, it took almost 23 hours to get the 
results, whereas with EXPOKIT it took less than 8 min- 
utes on the same machine. 



III. NUMERICAL RESULTS: COMPARISON 
BETWEEN INTEGRABLE AND 
NONINTEGRABLE MODELS 

We consider L even and choose as an initial state a 
sharp domain wall where the spins point up in the first 
half of the chain and down in the other half, 1^(0)) = 
| ftt • • ■ t-J- • • • ill) ■ This state belongs to the largest 
subspace, S z = 0, which has dimension dim = ( L y 2 ) = 
L!/[(L/2)!] 2 . 

Experimentally such a domain wall state may be cre- 
ated by the application of a strong spatially varying mag- 
netic field. The time-evolution is initiated when this 
strong field is quenched and the couplings become sud- 
denly effective. The system dynamics is then dictated by 
Hamiltonian (TTJ) with e„ = for all sites or, in the case 
of an impurity, with a small field remaining only on site 
L/2+1. 

The decision to place the impurity at site L/2 + 1 and 
not at L/2 has a reason. Placing it at L/2 implies the 
same results in terms of chaos indicators, but the trans- 
port of the excitations is impeded even before A surpasses 
1. Due to the Ising interaction and border effects, the do- 
main wall has a very large energy of(L — 3)A/4. As A 
increases, and especially when it becomes larger than 1, 
it becomes difficult to find states that are resonant with 
1^(0)) and the domain wall hardly moves. This localiza- 
tion occurs for even smaller values of A if the impurity 
is at L/2, because it adds to the already large energy of 
the domain wall state an additional positive contribution 
from the impurity field (23j . 

For the XX model, the width of the domain wall in- 
creases linearly in time |27j . As a result, the rightmost 
excitation (the up-spin in the middle of the domain wall) 
reaches the border of the chain at time t ~ L/ (2 J). Thus, 
to avoid border effects, we will focus on t < L/(2J) for 
models with NN couplings only, even in the presence of 
interaction and disorder. In the case of systems with both 
NN and NNN couplings, the excitations move faster (as 



can be seen in Fig. [3]) , so we will center our attention on 
times not much larger than i/(4J). 

A. Transverse spin correlation function 

We study the transverse spin correlation function, 
Eq. (|5|), between site L/2 and sites L/2 + m after a time 
long enough for the domain wall to have broadened and 
short enough to prevent effects from the borders. 

It was analytically shown for the XX model that the 
correlation function at t — > oo reaches a nonequilibrium 
steady state given by [27| , 

C xx (n, n + m) *±T C xx (m) cos (j^j , (10) 

where the ground state correlation function is, 

C xx {m) « ^Lk 2 , (11) 
V8m 

with \ogn = \ / °° f (e' 4t - ^r- t ). Thus the ground 
state correlation function is modulated by oscillations at 
a wavelength A where A = 2/po, ±Po bein g th e polariza- 
tion at the two halves of the domain wall [26|, [27| ■ From 
Eq. ([5]) and Eq. ©, these oscillations also imply a cur- 
rent flow of magnitude /iVi = {J /ft) sm(7rpo)- In our case 
Po=l/2. 

The equilibrium value of the xx-corrclation function, 
C xx (m) , is shown with a solid (black) line in all panels of 
Fig. [I] The absolute value of C xx (L/2, L/2 + m) for the 
XX model is depicted with circles in panels (a) and (b), 
and perfectly agrees with the analytical result. For the 
XXZ model, as A increases, the results deviate gradually 
from Eq. (fT0"|) . as illustrated in panel (a). For A > 0.1, it 
reaches a point where longer wavelengths become evident 
[panel (b)]. As will be shown in the next sub-section, 
this increase in the wavelength can be associated with a 
decrease in the current as A is increased. 

The results for the chaotic systems are shown in panels 
(c) and (d). The overall decay of the correlation func- 
tion is found to be more rapid than for the XX chain 
[Eq. (|10p]. Moreover, both chaotic models show an over- 
all increase in wavelength as A increases. As mentioned 
before, for the NN+NNN model, the simple correspon- 
dence between the nearest-neighbor correlation function 
(C xx (n,n + 1)) and current no longer holds. Thus even 
though the wavelength of oscillations increases for the 
NN+NNN model as A is increased from to 0.3, in this 
regime the current in fact also increases, as shown in 
Fig. 2 (c). 

B. Local spin current 

Spin and energy currents, I s and Ie respectively, have 
been extensively studied in the context of intcgrable pe- 
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FIG. 1: (Color online.) Absolute value of the transverse 
spin correlation function, \C xx (L/2, L/2 + m)\ vs m for time 
t = 10/J when only nearest-neighbor couplings are present 
and for t = 6/ J when next-nearest-neighbor couplings are 
also present; L — 20. For the three models, circle: A = 
and square: A = 0.3. Panels (a) and (b): integrable model 
with £_l/2+i = 0.0, a = 0.0. In panel (a) the results deviate 
monotonically from the case where A = 0.0 as the interaction 
increases: A = 0.02, 0.04, 0.06, 0.08, 0.10. Panel (c): chaotic 
model with £l/2+i = 0.5 J, ot — 0.0. Panel (d): chaotic model 
with £l/2+i = 0.0, a = 0.4. 



riodic chains [l|, 0, HH, [44[ ■ Ie commutes with the Hamil- 
tonian of the XXZ system implying divergent thermal 
conductivity. Although I s is conserved only for the XX 
model, it has been argued that spin transport should also 
be anomalous in the XXZ case when S z ^ [j], [ll[ • 

Numerical results for the local spin current, Eqs. © 
and (J7J, for the three models studied here are shown 
in Fig. [5] For NN couplings only, I s decreases as the 
strength of the Ising interaction increases for both the 
integrable model [panel (a)] and the chaotic case with 
impurity [panel (b)]. For the integrable model, this be- 
havior is expected. Since the energy of the domain wall, 
(L — 3) A/4, increases with A, the number of states capa- 
ble of coupling to it decreases and therefore the motion 
of the excitations becomes more restricted. Equivalcntly, 
the eigenstates become spatially more localized. For the 
system with impurity however, the addition of interaction 
leads to the onset of chaos and corresponds to eigenstates 
that are spatially more delocalized. (The relationship 
between the level of delocalization and the magnitude 
of the spin current is further explored in Sec. IIV"|) . At 
first sight, the existence of more delocalized eigenstates 
seems to contradict the decrease of the current, but this 
decrease may be due to a crossover from ballistic to dif- 
fusive transport. 

The transport behavior of the three models considered 
here have been studied before. It is ballistic in the inte- 
grable domain, which occurs when a = 0, A < 1, e„ = 
E S H and also whcn a = 0, A = 0, e L/2 +i ^ [Hj, 
but it becomes diffusive for the chaotic models. The re- 
lationship between non-ballistic transport and the onset 




FIG. 2: (Color online.) Local current, I S: „ [Eqs. (JBJ, (0)] vs 
n for time t — 10/J when only nearest-neighbor couplings 
are present and for t = 6/J when next-nearest-neighbor cou- 
plings are also present; L = 20. Panels (a): integrable model 
with e_L/2+i = 0.0, a = 0.0. Panel (b): chaotic model with 
£ l/2+i = 0.5J, a = 0.0. Panel (c): chaotic model with 
£l/2+i = 0.0, a = 0.4. Circle: A = 0; square: A = 0.5; 
up-triangle: A = 1.0; down-triangle: A = 2.0. 



of chaos in the single impurity model with interaction 
was investigated in Ref. [T(|. There it was shown that 
spin and thermal conductivity relax within a time that 
scales with the size of the chain. In the case of NNN 
couplings, the openingof a gap has been associated with 
diffusive transport 0,13, H Efl , but not much work exists 
that explicitly connects transport behavior with the level 
of chaoticity of the system. 

The model with NNN couplings and a — 0.4 is chaotic 
for any value of A < 1. The results for local current 
for this system, shown in panel (c) , are yet more surpris- 
ing. Here, the behavior of the current with A is non- 
monotonic. It first increases with A, reaching a maxi- 
mum level at A ~ 0.5, and only then starts diminishing 
with increasing A. The level of spatial delocalization, 
on the other hand, decreases monotonically with A (sec 
discussion in Sec. IV) and thus cannot justify the result 
for current. A more thorough analysis of the level of 
chaoticity in the presence and absence of Ising interac- 
tion, as well as how it relates to transport behavior, may 
shed some light on this scenario. 



C. Local magnetization 

The non-monotonic behavior of the model with a^0 
is reflected also in the results for the local magnetization 
in Figs. [3 H and El 

The contour plots in Fig. [3] show the absolute value of 
the magnetization of each site vs time. The interaction 
strength increases from top to bottom. For the integrable 
model (left panels), as A increases, the motion of the ex- 
citations becomes more difficult and it takes longer for 
the sites to reach very low absolute values of the mag- 
netization. A more complicated and asymmetric picture 
emerges in the presence of an impurity (middle panels), 
but the trend is still the same. However, when NN and 
NNN couplings are considered (right panels), as A in- 
creases from zero to 0.5, the excitations spread much 
more efficiently and lead to an enhanced decay of the 
absolute values of local magnetization. As A is increased 
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beyond 0.5, the trend finally reverses and the behavior 
becomes more similar to the other models, where the de- 
cay slows down with the Ising interaction. 
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FIG. 3: (Color online.) Contour plots of the absolute value of 
the magnetization of each site \M n \ [Eq. ([8])] vs time; L — 20. 
Left column: integrable model with e L / 2 +i = 0.0, a = 0.0. 
Middle column: chaotic model with cl/2+i = 0.5 J, a = 0.0. 
Right column: chaotic model with £l/2+i = 0.0, a = 0.4. 
Top panels A = 0; middle panels: A = 0.5; bottom panels: 
A = 1.0. The contour lines range from 0.05 in the darkest 
region to 0.5 in the lightest region. 
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FIG. 4: (Color online.) Magnetization of the first half of 
the chain [Eq. (J5J] vs time; L = 20. Panel (a): integrable 
model with cl/2+i = 0.0, a = 0.0. Panel (b): chaotic model 
with e_L/2+i = 0.5 J, a = 0.0. Panel (c): chaotic model with 
££/2+l = 0.0, a = 0.4. Solid line A = 0; dashed line: A = 0.5; 
dot-dashed line: A = 1.0. 

The local magnetization of the first half of the chain 
[Eq. ([9])] is shown in Fig. [4] and reinforces the observa- 
tions from the contour plots. The decay of M^t/i) is 
linear and the most rapid when A = in the case of the 
integrable clean model [panel (a)] and the model with a 



single impurity [panel (b)]. For these two cases, the de- 
cay slows down as the Ising interaction increases. Panel 
(c), on the other hand, provides a clear illustration of 
the non-monotonic behavior of the NN+NNN model with 
A. Notice the abrupt and linear decay of M( 1;i / 2 ) when 
A = 0.5, which is to be contrasted with the slow and 
nonlinear decay for A = and 1. 

Another perspective to the results of the contour plots 
is provided in Fig. [5l which shows the behavior of the 
magnetization M n of site n in time. From this figure one 
verifies the existence of two different velocities at play. 
One is related to how fast the magnetization can reach 
very small values. For the integrable model, M n reaches 
the lowest values faster when A=0, which justifies the 
steep slope of Mnr/2) (solid line) in Fig. [4] (a), whereas 
for the chaotic model with NNN couplings, this speed is 
higher for A=0.5. 

The other velocity is related to the motion of the do- 
main wall front. For the integrable case, this velocity 
increases with A. This is a well known result from Bethe 
Ansatz [46[ and bosonization 47 1 , and has also been stud- 
ied numerically in [48| . A rough idea of the front veloc- 
ity may be obtained by looking at the first moment the 
magnetization reaches its minimum value. For the in- 
tegrable model this happens earlier as the Ising interac- 
tion increases. Compare, for example, the time when site 
10 reaches the lowest magnetization for A = [top left 
panel] and A = 1 [bottom left panel]: it is after 2/ J for 
the first and before 2/ J for the latter. The same pattern 
is seen for all other sites. By passing an imaginary line 
through these minima, one sees that the slope becomes 
steeper as A increases. In the case of the chaotic model 
with NNN couplings, an evident feature is the suppres- 
sion of oscillations, but an estimate for the front velocity 
becomes more difficult. Focusing on the magnetization 
of the middle site 10, the velocity appears to also show 
a non-monotonic behavior. Contrary to the integrable 
case, the front velocity appears to be slightly larger for 
A = than for A = 0.5; it is only for A > 0.5 that the 
behavior for the two models become similar and the spin 
wave velocity increases again with interaction (49j . 



IV. DELOCALIZATION MEASURE 

In an attempt to understand the non-monotonic be- 
havior of the model with NN and NNN couplings, we 
study here the so-called inverse participation ratio (IPR) • 
For an eigenstate written in the basis vectors \<f>k) 



IPR is defined as 



IPR, EE 



E 



dim 
fe=l 



„fc|4 



(12) 



This quantity measures the number of basis vectors that 
contribute to each eigenstate, that is, it measures the 
level of derealization of each eigenstate HH . It is 
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FIG. 5: (Color online.) Magnetization of site n vs time; 
L — 20. Left panels: integrable clean model with ex/2+i = 
0.0, a = 0.0. Right panels: chaotic clean model with NNN 
couplings, e_L/2+i = 0.0, a = 0.4. In each panel, the curves 
range from site 10 to site 6, as indicated in the top left panel. 
Site 10 corresponds to the center of the domain wall; the up- 
spin placed here is the first one to move. Top panels A = 0; 
middle panels: A = 0.5; bottom panels: A = 1.0. 



small when the state is localized and large otherwise. No- 
tice that the result depends on the choice of basis. This 
choice is made according to the question one is after. In 
studies of spatial localization, for example, the site-basis 
is the most appropriate one. When the goal is to separate 
regular from chaotic behavior, the mean- field basis [5l| . 
corresponding to the eigenstates of the integrable limit, 
becomes more relevant. 

We consider here the site-basis, since we are interested 
in quantifying the spatial broadening of the domain wall. 
For the three models and different values of the Ising 
interaction strength, Table fl] gives the value of IPR av- 
eraged over all eigenstates, (IPR); the number of states 
that participate in the evolution of the initial domain 
wall state, = £) . Cj\if)j)e^ iE ^, with probability 

| Cj- 1 2 larger than 10 -6 ; and the IPR averaged over these 
latter eigenstates, (IPR)|c j |2 >10 -6. The purpose of these 
quantities is the following: (IPR) quantifies the overall 
level of derealization of the whole system taking into ac- 
count all L\/[(L/2)l] 2 eigenstates; the number of states 
with \Cj\ 2 > 10~ 6 identifies the eigenstates that most 
contribute to the evolution of the domain wall state; and 
finally (IPR)| Cj |2 >10 -6 indicates whether these contribut- 
ing eigenstates have or have not the same level of der- 
ealization of the whole system. 

For the three models, the behavior is the same as A 
increases from 0.5 to 1.0: the eigenstates become more 
localized. This is reflected in the dynamics of the three 
systems: the excitations become more localized, which 
decreases the magnitude of the spin current and slows 
down the decay of local magnetization (while the ve- 



TABLE I: Value of IPR averaged over all eigenstates; number 
of states participating in the evolution of the domain wall 
state, |*(t)) = J2j Cj\^j)e~ lEjt , with probability larger than 
10" 6 ; IPR averaged over these eigenstates that contribute 
with \Cj\ 2 > 10 -6 . System size: L = 16, subspace S z = 0. 



Model 
{a,e L/2+1 /J) 


A 


(IPR) 


|Ci| 2 's > 10" 6 


(IPR)|o,p> 10 -6 


(0.0,0.0) 


0.0 


3075.56 


3054 


3009.98 


(0.0,0.0) 


0.5 


2358.02 


2211 


2333.44 


(0.0,0.0) 


1.0 


1732.02 


1139 


1430.35 


(0.0,0.5) 


0.0 


1921.53 


6696 


1917.91 


(0.0,0.5) 


0.5 


3679.49 


6222 


3737.60 


(0.0,0.5) 


1.0 


3118.76 


2868 


2439.36 


(0.4,0.0) 


0.0 


3906.90 


4378 


4011.37 


(0.4,0.0) 


0.5 


3684.89 


2241 


3624.25 


(0.4,0.0) 


1.0 


2840.99 


663 


1670.35 



locity of the domain wall front increases). However, as 
A changes from 0.0 to 0.5, this correspondence remains 
valid only for the integrable system, whereas for the other 
two models orthogonal behaviors are seen: 

(i) In the case of the system with impurity, the ad- 
dition of interaction is followed by a crossover to chaos, 
which justifies the substantial increase of the value of 
IPR. The fact that the magnitude of the spin current 
decreases, going contrary to derealization, may be indi- 
cating a transition from ballistic to diffusive transport, 
which should slow down the motion of the excitations. 

(ii) The system with NN and NNN couplings is chaotic 
when A = 0.0 and also when A = 0.5. Based on the 
spatial contraction of the eigenstates that occurs when 
the Ising interaction is turned on, one might expect a 
decline in the motion of the excitations. However, what 
is verified in Figs. [21 El El and [5] is just the opposite. The 
interpretation of this counter-intuitive behavior requires 
further studies. 

It is interesting to compare the values of (IPR) and 
(IPR)|C |2>io- 6 - They are very similar for small A, but 
differ when A = 1. As mentioned before, the energy of 
the domain wall increases with the strength of the Ising 
interaction. At A = 1 the energy is already very large 
and close to the edge of the spectrum. In systems with 
few-body interactions (ours have only two-body interac- 
tions) , eigenstates with energ y cl ose to the edges of the 
spectrum are more localized [52l [53| . Thus, the eigen- 
states that contribute to the evolution of the domain wall 
state are fewer (the numbers in the column of \Cj\ 2 > 
10~ 6 decrease significantly from A = 0.5 to A = 1.0) and 
also more localized than those participating in the evolu- 
tion of generic initial states ((IPR)|c , j |2 >10 -6 < (IPR)). 
This should affect the relaxation process of the domain 
wall state, which should become slower for large A. 
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V. CONCLUSION 

We considered a one-dimensional spin- 1/2 system in 
both integrable and chaotic regimes. The intcgrablc limit 
corresponded to a clean chain in the presence of nearest- 
neighbor couplings only. Chaos was induced (i) by adding 
a single-site impurity field or (ii) by adding next-nearest 
neighbor couplings. The time-evolution of an initial state 
in the form of a sharp domain wall was studied. 

The analysis of the dynamics focused on transverse 
spin correlation function, local spin current, and local 
magnetization. In order to understand the behavior of 
these quantities, we resorted to derealization measures; 
the motivation being that usually a good picture of the 
static properties of a system may indicate what to ex- 
pect for its time evolution. This approach served us well 
when dealing with the integrable model. Here, as the 
Ising interaction strength A increases, the level of spatial 
derealization of the eigenstates lowers and, as one might 
expect, the magnitude of the spin current decreases and 
the decay of the absolute value of local magnetization 
slows down. 



The same behavior was verified for the chaotic models 
for the case of large A, but for small Ising interaction, 
the level of derealization changed in a way contrary to 
the change in current and magnetization. For model (i), 
the addition of interaction leads to the onset of chaos 
and consequently the eigenstates further dclocalize. The 
magnitude of the current, however, decreased. We spec- 
ulate that this is due to the transition from ballistic to 
diffusive transport. For model (ii), the system is already 
chaotic at A = and the eigenstates only shrink with 
the Ising interaction. Contrary to that, the motion of 
the excitations is enhanced, as attested by the results for 
spin current and local magnetization. Further studies are 
needed to clarify this behavior. 
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